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We have studied electron transport in clean single-walled carbon nanotube quantum dots. Because 
of the large number of Coulomb blockade diamonds simultaneously showing both shell structure and 
Kondo effect, we are able to perform a detailed analysis of tunneling renormalization effects. Thus 
determining the environment induced level shifts of this artificial atom. In shells where only one 
of the two orbitals is coupled strongly, we observe a marked asymmetric gate-dependence of the 
inelastic cotunneling lines together with a systematic gate dependence of the size (and shape) of 
the Coulomb diamonds. These effects are all given a simple explanation in terms of second-order 
perturbation theory in the tunnel coupling. 



A solid-state quantum dot (QD) is often referred to as 
an artificial atom[l|, and the energy levels of this ar- 
tificial atom can be measured usingbias spectroscopy in 
the Coulomb blockade regime 0, Q. In this regime, 
the couplings to the metallic leads are usually regarded 
as small perturbations, giving rise to a finite conductance 
without otherwise affecting the internal structure and en- 
ergy levels of the QD. Deviations from the perfect shell 
structure in, for example, carbon nanotubes are usually 
attributed to imperfections of the single-walled carbon 
nanotube (SWCNT). Consequences of stronger tunnel 
coupling to the metallic leads have been studied in the 
context of broadening of energy levels Q and Kondo ef- 
fects 

For moderate coupling strengths a certain degree of 
level shifts and broadening is expected to arise from the 
hybridization with the leads. Nevertheless, these effects 
are difficult to distinguish and therefore rarely discussed 
in the literature. So far, the only case where such renor- 
malization effects have been shown to have important 
observable consequences is for QD spin-valve systems, 
where either non-coUinear magnetization gives rise to 
off- diag onal renormalization of the spin structure on the 
QD 111 or spin-dependent tunneling gives rise to a gate- 



1^ , isj . In a similar 



dependent exchange-field on the QD 
manner, one might expect a difference in couplings to the 
two orbitals (or subbands) in the SWCNT to give rise to 
distinct observable renormalization effects in a SWCNT. 

In this Letter, we report on measurements on sev- 
eral SWCNT-devices showing regular four-electron shell 
structure. The best device shows a total of 60 regular 
four-electron shells, allowing for a unique quantitative 
and statistical analysis of renormalization effects. In par- 
ticular, we study the cases where the four-electron shells 
exhibit Kondo effect in only one of the diamonds in the 
quartet, indicating that one of the (almost) degenerate 
orbitals has a stronger coupling to the leads than the 
other. Such behavior is seen in 28 four-electron shells. 
Based on this analysis we find that the relative mag- 
nitudes of the tunnel coupling of orbitals one and two 
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FIG. 1: Bias spectroscopy plot of a SWCNT QD for -lOV < 
Vg < lOV. Coulomb diamonds are seen for almost every 
added electron (285), and 88 odd-occupancy Coulomb block- 
ade diamonds exhibit a zero-bias Kondo resonance. 



determine three distinctly different types of quartets. 

A difference in tunnel coupling for the two orbitals 
gives rise to two features: 1) the inelastic cotunneling 
threshold given by the orbital splitting develops a marked 
gate dependence, and 2) the addition energies of the indi- 
vidual Coulomb diamonds within a quartet deviate from 
the constant-interaction model , depending on the rela- 
tive magnitude of the tunneling couplings. Both features 
can be explained by level shifts of the many-body states 
of the QD, calculated here to second-order in the tunnel 
coupling. 

The devices were made on a highly n-doped 500/im 
thick silicon substrate with a 500nm thick Si02 top layer. 
Chemical vapor deposition (CVD) was used to grow the 
SWCNTs from predefined catalyst particle areas. The 
SWCNTs were contacted by thermal evaporation of tita- 
nium/ aluminum/ titanium (5nm/ 40nm/ 5nm) plj . and 
the electrodes were separated by 300nm and patterned 
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FIG. 2: Differential conductance as a function of bias voltage, Vad, and gate voltage, Vg (left). Schematic of the four-electron 
shells (right), (a) Four-electron shell filling for Fi ~ with strong enough coupling to facilitate kondo ridges in the N+1 and 
N+3 diamonds, as shown in the schematic. A total of 17 four-electron shells show this type of behavior across the full gate 
voltage range, (b) Four-electron shell filling for Fi <^ F2. In the N+1 and N+2 diamonds, the inelastic cotunneling threshold 
has clearly acquired a gate-dependence. In total there are 21 four-electron shells showing this type of behavior. Here the 
right panel shows the cotunneling-thresholds and diamond edges as calculated within second-order FT (The white Kondo-ridge 
is drawn by hand), (c) Four-electron shell filling for Fi S> F2. Gate-dependent cotunneling thresholds are now observed in 
the N+2 and N+3 diamonds instead. In total there are 7 four-electron shells showing this type of behavior. 15 four-electron 
shells could not be uniquely categorized. In the 28 shell-sequences where Fi 3> F2 or Fi <C F2 the gate-dependent cotunneling 
threshold ridges inside the diamonds appear to form a U or V lying down with the openings facing the zero-bias Kondo ridge. 
Together the two U's (V's) and the Kondo ridge resemble a two-headed arrow, where the U's (V's) form the two heads and the 
Kondo ridge the shaft. This means that the tree cente r diamon ds in the di fferent d omains have co tunnelin g and Kondo ridges, 
which can be schematically summarized as: Fi ~ r2 : — = — , Fi ^ F2 : << — and Fi S> F2 : — >> 



by electron beam lithography. The silicon substrate was 
used back-gate to change the electrostatic poten- 
tial of the SWCNT QD. The electron transport measure- 
ments were carried out in a cryogenic insert system with a 
base temperature of 300mK. Two-terminal conductance 
was measured using lock-in techniques with 13^V AC ex- 
citation amplitude and a DC bias voltage (Vgd) applied 
to the source, while the drain was grounded through a 
current amplifier. 

Room temperature measurements of the conductance 
as a function of back gate voltage {Vg) indicate that the 
SWCNT in the presented device is metallic. In the full 
gate voltage span, — lOV <Vg< lOV, four-electron shell- 
filling is observed at low temperatures, and an overview 
of the regular behavior of the device can be seen in Fig. [1] 
In total, 60 full four-electron shells are observed across 
the measured gate voltage range. Disregarding a few 
gate-switches, we have observed 285 consecutive electron 
additions. 60 four-electron shells can be identified, the 
remaining 45 electrons are counted in incomplete shells. 



We have extracted the Kondo temperatures from 88 equi- 
librium Kondo resonances from the width of the Kondo 
peak {Wfwhm ^ 4A:b7V)[1J|. The Kondo tempera- 



tures found by measuring the Kondo peak widths arc in 
good agreement with Kondo temperatures found for se- 
lected peaks, when measuring the zero-bias conductance 
as a function of temperature, and fitting to the universal 
scaling function ^15l| (not shown). We observe no system- 
atic gate-dependence of the Kondo temperatures, with 
an average of (Tk) ~ 1.2K and standard deviation 0.5K. 

A more detailed bias spectroscopy plot of one of the 
four-electron shells from Fig. [1] can be seen in Fig. [2] (a). 
The plot resembles previously reported characteristics of 
a SWCNT QD in the Kondo regime [i,^. In diamonds 
N+1 and N+3, corresponding to odd electron occupa- 
tion, spin-1/2 Kondo ridges can be seen as horizontal 
lines at zero bias voltage. In the N+1, N+2, and N+3 
diamonds, the onset of inelastic cotunneling can be seen 
as practically horizontal lines at Vsd ~ ±0.7mV [lit . 
The cotunneling ridge at Vsd ~ +0.7mV in the N+2 dia- 
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FIG. 3: (Color online) The ratio of the widths of the dia- 
monds (proportional to the addition energies) in the 28 un- 
usual four-electron shells to the odd-electron diamonds with- 
out a Kondo ridge in the same four-electron shells. The solid 
(cross-hathed) pairs is for Fi <C F2 (Fi S> F2). The average 
width of the even-electron diamonds (black) is 1.47, and the 
average width of the odd-electron diamonds with a Kondo 
ridge (red (dark grey)) is 0.66, as indicated by the black lines. 

mond is significantly Kondo enhanced [lo| . Equilibrium 
and nonequilibrium Kondo features are observed in the 
full gate voltage span. From the plot we can extract 
the characteristic energies of the system as indicated in 
the schematic drawing on the right in Fig. [2] (a) 0, \v\ . 
We obtain a Coulomb energy oi U = 2.9meV, orbital 
splitting 5 — 0.7meV and single-particle level spacing 
A£; = 3.5meV. 

In addition to this well-established four-electron shell 
behavior, we have observed two distinctly different 
types of four-electron shells, exemplified by the plots in 
Figs. [21(b) and (c). In Fig. [21(b), there is a Kondo ridge 
in N+3, but the expected Kondo ridge in the N-l-1 di- 
amond is absent. In addition, the cotunneling thresh- 
olds in the N-t-1, and N-f2 diamonds increase from left 
{Vsd ~ ±QAmV) to right (T^d ~ ±1.2toF), showing al- 
most identical gate-dependence. In Fig. [2l (c), the sit- 
uation is reversed: A zero-bias Kondo ridge now only 
occurs in the N-|-l diamond, and the cotunneling thresh- 
olds in the N-l-2, and N-l-3 diamonds decrease from left 
{Vsd « ±QAmV) to right [Vsd ~ OmV). The electron- 
hole symmetry (around the middle of the N+2 dia- 
mond Iql) which is observed in panel (a) is clearly lifted 
in panels (b) and (c). Furthermore the main features in 
panels (b) and (c) are each others mirror image. 

As already mentioned, the sample shows a large num- 
ber of four-electron shells (60). 28 of them show only 
a single equilibrium Kondo resonance and exhibit gate 
dependent inelastic cotunneling lines inside the two ad- 
jacent Coulomb diamonds. 

Furthermore, there is a clear correlation between the 
size of the odd-electron diamonds with, and without 
the zero-bias Kondo ridge, respectively. This is demon- 



strated in Fig. [3l where the widths of the A'' + 2 dia- 
monds (black) and the odd-electron diamonds exhibiting 
a Kondo resonance (red (dark grey)) are shown relative 
to the widths of the odd-electron diamonds not showing a 
Kondo resonance. For all 28 quartets of this type, the di- 
amond showing a zero-bias Kondo ridge is narrower than 
the diamond without the Kondo ridge, independent of 
whether Fi < r2 or Fi > F2. Observing this in 28 out 
of 60 shell-sequences, this behavior cannot be dismissed 
as a measurement irregularity. We have seen similar gate 
dependences in several other SWCNT devices with inter- 
mediate coupling to the electrodes. 

We shall now focus on the quartet mapped out in 
Fig. [21(b), where Fi <C The presence and absence of 
a zero-bias Kondo-ridge in the N+3, and N+1 diamonds, 
respectively, already indicates that the two orbitals must 
have different tunnel couplings, since the Kondo tempera- 
tures must be respectively larger than or smaller than the 
base temperature (300 mK). The Kondo temperatures 
are given by T^!^ - V4rWexp(~7r[//(8F,) [ill, in terms 
of the individual tunnel-broadenings F,; = '^a=s d^i,on 
with Vi a = TTi^Fl^i.Qp, determined by the density of 
states in the source, and drain electrodes, vp, and their 
respective tunnel couplings to the dot, ti^a, assumed to 
be energy independent in the relevant energy range. The 
lack of particle-hole symmetry in this quartet is thus con- 
sistent with the orbital symmetry being lifted by a dif- 
ference between Fi and F2. The F^ are generally diffi- 
cult to extract, but by fitting the even valley part of the 
Coulomb peaks of Fig. [21(b) to a Lorentzian (cf. Ref . fiol) . 
we estimate that Fi « 0.08mcV and F2 « 0.46meV. This 
yields T^-* ~ 7/iK and T^-* ~ 2K, consistent with the ob- 
servation of a Kondo ridge only in the N+3 diamond at 
T - 300 mK. 

The finite-bias cotunneling lines inside the N+1 and 
N+2 diamonds of Fig. [2l (b) can be ascribed to an inelas- 
tic cotunneling processes in which an electron traversing 
the dot from source to drain has enough energy to excite 
a single electron from orbital one to orbital two. The in- 
dividual many-body states of the dot electrons are renor- 
malized by virtual charge-fluctuations, and, as we shall 
argue below, these tunncling-induced level shifts can in- 
deed give rise to a gate-dependent cotunneling threshold. 

Within a single four-electron shell, the many-body 
states of the dot electrons can be enumerated by the 
number of electrons occupying the two highest lying 
single-particle orbitals. In terms of bare level-position, 
Ed, intra-dot Coulomb repulsion, [/, and orbital split- 
ting, 5, the corresponding energies are given by Ei^j = 
{i + i)£d + i5 + j - l)C//2, where i{j) de- 

notes the number of electrons in orbital 1(2). The level 
shifts, SEij = Eij — Eij, are determined within second- 
order many-body perturbation theory (FT) in the tunnel 
coupling, by considering all possible fluctuations experi- 
enced by a given charge conflguration [20| . For example. 
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the charge state |1, 0) connects to |0, 0) and to |2, 0) and 
|1, 1), respectively, by virtual tunneling out^ and in of an 
electron. This leads to the following shift in energy: 



CK — S.d 



— Re 

D 



TT 



+ 



ria(l - fjiO - fJ-a)) 



(1) 



Ein + UJ — i?2,o + ir Ei,o + uj - EiA + iV 



where 2D is the conduction-electron bandwidth, f{Lo) 
is the Fermi function, and )jLs,d the chemical potentials 
of the source, and drain electrodes. A level broadening 
of order F ~ Fi + r2 has been included in the energy 
denominators and serves to cut off the otherwise loga- 
rithmically singular result. Notice that the last process 
has multiplicity two from spin, and that only this pro- 
cess has the amplitude F2. The excited state is renor- 
malized correspondingly and <5i?o,i is given by Eq. ^ 
with all orbital indices reversed, i.e. Ei,j — *■ Ej^i and 
Fia ^ F2Q. This means that for gate- voltages close to 
the 0/1 (left diamond corner), or 1/2 (right diamond 
corner), charge-degeneracy point, respectively, it is the 
first, or the two last terms in Eq. ([1]) which dominate the 
renormalization, with small energy denominators. When 
Fi <C F2, the i?o,i energy level (orbital 2) will have a 
large shift at the left corner of the diamond (0/1) and 
the Eifi energy level (orbital 1) will have a large shift 
at the right corner (1/2), as illustrated in Fig. [H This 
means that the cotunneling threshold increases from left 
to right when Fi ^ F2, and decreases from left to right 
when Fi ^ F2. To be more precise, the inelastic cotun- 
neling threshold inside the N-l-1 diamond is determined 
as the energy, 5, solving the equation 5 = Ei^ — i?04, 
with fXs — fj,d = ±6/e in the expressions for the Eij. The 
solution, with parameters found above [l^l, is plotted in 
the right panel of Fig. [2] (b), with a reasonable resem- 
blance to the gate-dependence observed in the left panel. 
Expanding near the middle of the A^ + 1 diamond, we find 
that dS/dsd « -4(F2-Fi)/(7rt7), for Sd = -U/2. Thus if 
F2 > Fi the cotunneling threshold clearly increases with 



gate voltage (eVg 



-^Ed), i.e. in the direction towards 



the stronger coupled -I- 3 diamond, as we observe. 

Calculating all the relevant level shifts, one can also 
determine the changed gate voltage width, AVg{n), of 
the N + n diamond. We find that AVg(l) « [/ - 
(2Fi/7r)ln(t7/F), AVg{2) ^ U + S and AFg(3) « [/ - 
(2F2/7r) ln(J7/F), strictly valid only in the regime where 
S ^ T ^ U . A larger S implies more correction terms, 
but the main trend remains the same: a large A^ -I- 2 
diamond and a small A^ -I- 3 diamond, as compared to 
the A^ 4- 1 diamond (Fi <C F2) in agreement with the 
statistics presented in Fig. [3l 

Finally, we have calculated the shifted Coulomb di- 
amond boundaries, using the second-order level shifts. 
For the found parameters, this gives rise to the diamond 




FIG. 4; (Color online) Schematic drawing illustrating the 
difference in renormalization of Ei^ (red (dark grey) solid 
lines) and _Bo,i (red (dark gray) dashed lines), at opposite 
ends of the A*" + 1-diamond. Assuming Pi <^ T2, the (0,1)- 
level is strongly perturbed by virtual emptying (left), whereas 
the (l,0)-level is most strongly perturbed by virtual double- 
occupancy (right). 



boundaries seen in right panel of Fig. [2] (b). The de- 
tails of this figure are not entirely reliable, and more 
work is needed in order to proceed beyond second-order 
FT and to incorporate the nonequilibrium mixing of 
ground, and excited states within each diamond. Nev- 
ertheless, this simple analysis indicates that tunncling- 
induced level shifts may give rise to skewed (i.e. non- 
parallel edges), rather than merely shifted and rescaled 
Coulomb diamonds, which is in fact seen in several of our 
experimental quartets, e.g. Fig. [21(c). 

In conclusion, we have fabricated gated SWCNT de- 
vices with intermediate transparency contacts. Trans- 
port measurements have revealed three distinct types of 
four-electron shells, two of which show a marked gate- 
dependence of the inelastic cotunneling ridges, together 
with a systematic pattern for the individual diamond- 
widths within a shell. Calculating second-order level 
shifts for the many-body states on the QD, we have ex- 
plained all of these observations in terms of a difference 
in tunnel coupling to the two orbitals in the SWCNT. 

We acknowledge support from EU projects ULTRA- 
ID, SECOQC, CARDEQ, and the Danish Research 
Council. 
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